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1.  Introduction 


Boron  suboxide  (BeO)  is  a  promising  material  system  for  lightweight  armor 
applications.  It  has  an  ultra-high  hardness1,2  (24-45  GPa),  low  density  (-2.61 
g/cc),  high  mechanical  strength  (26-30  GPa),3  high  oxidation  resistance4  (below 
1,200  °C),  and  is  chemically  inert.  However,  it  has  yet  to  be  fielded  in  Soldier- 
mounted  protection  due  to  several  processing  and  perfonnance-based  challenges 
including:  1)  trans-granular  failure  resulting  in  low  fracture  toughness,5,6  2)  a 
reduction  in  shear  strength  under  extreme  environmental  conditions7,8  such  as 
high  temperature  and  pressure,  and  3)  manufacturing  challenges  (i.e.,  narrow 
processing  window,  high  pressure)  required  to  create  highly  crystalline,  fully 
dense  stoichiometric  BeO.9  To  maximize  the  performance  of  BeO  as  a  lightweight 
armor  material,  we  must  first  understand  the  competing  influences  of  crystal 
structure,  defect  populations  (i.e.,  vacancies,  twins,  stacking  faults,  carbon  [C] 
contamination),  stoichiometry,  and  surface  reconstructions  on  the  nature  of  the 
atomic  bonding  within  the  bulk  crystal  and  near  interfaces  (i.e.,  grain/phase 
boundaries,  free  surfaces)  and  associate  these  effects  with  specific  failure 
mechanisms  (i.e.,  grain  boundary  decohesion,  twinning/stacking  fault  formation, 
intragranular  fracture,  and  bulk/grain  boundary  amorphization). 

B60  is  particularly  sensitive  to  carbon  contamination,  which  can  be  introduced 
during  any  ceramic  processing  stage  (for  example,  residual  organics  during 
powder  synthesis  due  to  precursors  and  ceramic  binders,  carbon-based  mixing 
media,  CO-containing  atmosphere  during  calcination  and  sintering,  graphite 
sintering  sleeves  and  furnace  components,  and  mechanical  treatments  such  as 
machining  and  polishing).  Our  previous  computational  study10  has  shown  that 
carbon  is  highly  stable  in  BeO.  The  precipitation  of  B4C-like  structure  within  the 
B60  unit  cell  is  exothennic,  and  the  energy  released  increases  as  more  carbon 
substitutes  into  the  lattice.  However,  to  date  a  complete  understanding  of  the 
influence  of  C  contamination  on  the  electronic  structure  and  mechanical 
properties  of  BeO  is  lacking.  In  this  report,  we  explore  these  relationships  by 
studying  the  effect  of  stoichiometry  and  C  contamination  on  the  nature  of  atomic 
bonding  (via  electron  localization  function),  the  electronic  structure,  and  the 
mechanical  properties  (i.e.,  elastic  stiffness  tensor,  bulk  modulus,  Young 
modulus,  and  shear  modulus)  of  bulk  BeO. 

2.  Computational  Details 

First  principles  density  functional  theory  calculations  using  the  generalized 
gradient  approximation  (GGA)  given  by  Perdew-Bruke-Emzerhof1 1  and  the 
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projector  augmented  wave  (PAW)  method12-14  were  performed  within  the  Vienna 
Ab  initio  Simulation  Package  (VASP).1'  The  plane  wave  cut-off  energy  was  set  to 
500  eV.  The  calculations  were  considered  converged  when  the  maximum  force  on 
each  ion  was  less  than  0.001  eV/A,  and  the  change  in  total  energy  was  less  than 
10-4  eV.  The  simulation  cell  consisted  of  a  42-atom  BeO  unit  cell  in  hexagonal 
representation  that  was  sampled  with  a  2  x  2  x  1  K-points  Monkhorst-Pack  grid  to 
obtain  structural  properties.  We  performed  a  full  relaxation  of  the  volume,  lattice 
parameters,  and  all  internal  atomic  coordinates.  The  calculated  equilibrium  lattice 
parameters  for  BeO  unit  cell  are  a  =  5.393  A  and  c  =  12.327  A,  which  are  in  good 
agreement  with  the  experimental  values  of  Higashi16  (a  =  5.374  A  and  c  =  12.331 
A)  and  Hubert17  (a  =  5.3902  A  and  c  =  12.3125  A,  c/a  ratio  =  2.284).  Hubert’s 
B60  was  oxygen  (O)  deficient  (maximum  achieved  O  occupancy  =  0.95  at  6 
GPa).  Therefore,  when  extrapolated  to  ideal  stoichiometry,  Hubert’s  lattice 
parameters  are  a  =  5.399  A  and  c  =  12.306  A.  Calculations  that  considered  the 
effect  of  carbon  doping  required  a  larger  2x2x1  supercell  (168  atoms)  to 
prevent  the  artificial  interaction  between  defect  states  across  the  periodic 
boundary.  The  lattice  parameters  calculated  at  the  T-point  (a  =  5.396  A  and  c  = 
12.31 1  A)  are  similar  to  the  3  x  3  x  3  k-points  calculations  (a  =  5.394  A  and  c  = 
12.3 18  A).  Therefore,  T-point  calculations  were  used  to  study  the  doped  system. 

3.  Results  and  Discussion 


3.1  Structural  and  Electronic  Properties 

To  understand  how  carbon  contamination  affects  the  bonding  character  and 
electronic  structure  of  BeO,  we  performed  a  detailed  analysis  of  the  total  density 
of  states  (DOS),  charge  density  distribution,  and  corresponding  electron 
localization  function  (ELF).  The  ELF  is  defined  by  Becke  and  Edgecombe18  as 
the  probability  density  for  finding  a  second  like-spin  electron  near  a  reference 
point.  ELF  is  a  dimensionless  relative  measure  of  the  electron  localization,  which 
is  calibrated  with  respect  to  a  probability  density  for  a  uniform  electron  gas.  High 
ELF  values  (ELF  ~  1)  indicate  that  the  electrons  are  more  localized  than  in  a 
uniform  electron  gas  of  the  same  density.  It  is  possible  to  differentiate  between 
ionic,  covalent,  and  metallic  bonding  from  a  qualitative  analysis  of  the  ELF  map. 
Covalent  bonding  is  indicated  by  a  maxima  in  the  ELF  occurring  along  the  bond 
midway  between  the  2  atoms.  Ionic  bonding  is  indicated  by  spherical  ELF 
distributions  around  the  atom  cores  with  area  along  bond  length  of  zero 
probability.  As  the  degree  of  covalency  increases  within  the  ionic  bond,  the 
attractor  (maxima)  migrates  towards  the  center  of  the  bond  and  does  not 
completely  circumscribe  the  atomic  cores.  The  relative  position  of  the  attractor 
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indicates  the  polarity  of  the  chemical  bond.  Metallic-like  bonding  is  indicated 
when  the  ELF  is  delocalized  over  a  region  that  is  not  covered  by  the  cores  (i.e., 
diffusely  distributed).  Figure  1  shows  the  total  charge  density  and  ELF  for  the 
undoped  BeO. 


Fig.  1  Total  charge  density  for  ILO  at  different  isovalues:  a)  1.2  (15%  of  max.),  b)  0.81 
(10%  of  max.),  c)  0.81  (another  view),  and  d)  0.54  (7%  of  max.).  ELF  at  different  isovalues: 
e)  0.89  (90%  of  max.),  f)  0.7  (70%  of  max.),  g)  0.7  (another  view),  and  h)  0.25  (25%  of  max.). 

As  discussed  in  our  previous  report,10  there  are  several  distinct  bond  types  within 
the  B60  structure:  1)  intra-icosahedral  bonds  (1.78  A)  between  boron  forming  the 
equator  of  the  icosahedra,  2)  intra-icosahedral  bonds  (1.76  A)  between  polar 
boron  at  the  icosahedra  apex,  3)  intra-icosahedral  bonds  (1.78-1.81  A)  between  a 
polar  and  equatorial  boron  within  the  icosahedra,  4)  inter-icosahedral  bonds 
(1.7  A)  between  2  polar  boron  atoms  from  different  icosahedra,  and  5)  inter- 
icosahedral  bonds  (1.5  A)  between  equatorial  boron  and  oxygen.  The  O  atoms 
occupy  the  interstitial  volume  between  the  icosahedra.  Each  O  atom  coordinates 
to  3  equatorial  Be  atoms,  each  from  different  neighboring  icosahedra. 

From  the  charge  density,  it  is  clear  that  the  bonding  between  boron  and  oxygen 
atoms  is  primarily  ionic.  As  illustrated  in  Fig.  la,  the  charge  density  isosurface  at 
isodensity  value  equal  to  1 .2  is  localized  around  the  oxygen  ions.  The  Bp  -Bp 
(subscripts  p,  e,  and  I  are  used  to  represent  polar,  equatorial,  and  interstitial  sites 
respectively)  inter-icosahedral  bonding  between  boron  atoms  connecting  2  B12 
icosahedra  are  shorter  and  have  stronger  covalent  character  than  the  bonds 
connecting  boron  within  the  icosahedra.  As  the  charge  density  isovalue  is 
decreased  to  10%  of  maximum  (Fig.  lb),  the  charge  density  lobes  show  highly 
directional  covalent  bonding  between  the  Bp  -Bp  bonds  joining  the  2  icosahedra. 
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Also  at  this  isovalue,  an  isosurface  ring  joining  the  6  equatorial  Be  atoms  begins 
to  develop  (Fig.  lc).  The  boron  atoms  within  the  B12  icosahedra  are  connected  to 
each  other  by  3  center  covalent  bonds.  As  the  charge  density  isovalue  decreased 
further  to  7%  of  the  maximum,  the  isosurface  spreads  to  the  entire  region  around 
the  B12  icosahedra  creating  a  “metallic  like”  surface  (Fig.  id). 

Figure  1  also  includes  the  ELF  for  this  system  at  various  isovalues.  The  ELF 
value  near  0.89  (Fig.  le)  highlights  the  electron  localization  along  the  Bp-Bp 
bonds  between  different  icosahedra.  As  already  noted,  this  bond  is  shorter  and 
consists  primarily  of  a  highly  directional  covalent  bond.  When  the  ELF  value 
decreases  to  0.7  (Fig.  If,  lg)  the  ELF  surface  encompasses  all  of  the  atoms.  This 
ELF  surface  shows  3  distinct  bonding  regions  (basins):  (Region  1)  shorter  Bp-Bp 
bonds  connecting  2  icosahedra,  (Region  2)  surrounding  O  atoms,  and  (Region  3) 
surrounding  the  B12  icosahedra.  Each  region  represents  a  distinct  bonding  type. 
Region  1  contains  the  stronger  highly  directional  covalent  bond  between  2  Bp 
atoms  in  2  icosahedra.  Region  2  includes  an  ionic  bond  between  an  equatorial 
boron  and  an  oxygen  atom  whereas  region  3  is  the  intra-atom  bonding  of  boron 
(Bp,  Be)  atoms,  which  make  up  the  icosahedral  frame.  As  the  ELF  value  is  further 
decreased,  the  3  basins  merge  into  each  other  (Fig.  lh).  We  performed  the  Bader19 
charge  analysis  for  this  system.  The  charge  transfer  from  polar  Bp  atoms  is  in  the 
range  of  0.01-0.09  electrons,  emphasizing  the  covalent  nature  of  the  bonds 
connected  with  these  atoms.  For  the  equatorial  Be  atoms,  the  charge  transfer  is  in 
the  range  of  0.45-0.48  electrons.  The  charge  gained  by  the  oxygen  atoms  is 
between  1.54  and  1.55  electrons.  Each  O  atom  is  bonded  to  3  Be  atoms  from 
different  icosahedra.  The  charge  transfer  to  the  O  atom  confirms  the  ionic  nature 
of  the  Be-0  bond. 

It  is  interesting  to  note  the  behavior  of  the  charge  density  and  ELF  after  carbon 
doping  (see  Appendix  A).  When  C  is  doped  at  an  oxygen  site,  the  sp2  hybridized 
carbon  forms  bonds  with  3  neighboring  equatorial  Be  atoms  (Be-C  bond  length  is 
1.52  A).  The  resulting  B-C-B  bond  angles  are  119.39°,  which  is  nearly  equal  to 
the  standard  120°  bond  angle  for  sp2  hybridization.  The  charge  density  isosurface 
is  evenly  distributed  around  the  C-Be  bonds,  whereas  the  ELF  shows  3  lobes 
along  3  C-  Be  bonds  indicating  strong  covalent  bonding  between  these  atoms.  If  2 
carbon  atoms  are  doped  at  2  neighboring  oxygen  sites,  both  carbon  atoms  become 
sp2  hybridized  and  fonn  bonds  with  nearest  neighboring  equatorial  boron  atoms. 
When  an  interstitial  is  placed  between  these  2  carbon  atoms,  forming  C-O-C,  C- 
B-C  or  C-C-C  chains,  the  bonding  changes  significantly.  If  the  interstitial  is  either 
carbon  or  boron,  the  2  tenninal  carbon  atoms  become  tetrahedrally  “sp3  like” 
bonded,  whereas  the  interstitial  boron  or  carbon  atom  are  coordinated  only  to  the 
tenninal  carbon. 
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The  variation  in  the  bond  angles  around  the  carbon  atom  reflects  the  changing 
nature  of  the  bonding.  For  the  C-B-C  chain  configuration  (Co-Be  and  Co-Bi  bond 
lengths  are  1.57  A  and  1.44  A,  respectively),  the  3  Be-C-Be  bond  angles  are 
116.6°,  whereas  the  Be-C-Bi  bond  angles  are  100.75°.  For  the  C-C-C  chain 
configuration  (Co-Ci  and  Be-Co  bonds  distances  are  1.32  A  and  1.61  A, 
respectively),  the  3  Be-Co-Be  bond  angles  are  1 14.49°,  whereas  the  Be-Co-Cibond 
angles  are  103.82°.  This  indicates  that  the  C-C  bond  is  stronger  than  the  Be-C 
bond.  Also,  the  Be-C  bonds  weaken  for  C-C-C  chain  when  compared  with  the  C- 
B-C  chain.  This  is  reflected  in  the  larger  Be-C  bond  distances  in  C-C-C  chain 
(1.61  A)  when  compared  with  the  C-B-C  chain  (1.57  A).  Both  the  B-C-B  and  C- 
C-C  chain  remain  linear.  However,  the  C-O-C  chain  is  slightly  bent  with  the  angle 
150.14°.  The  other  Be-C-Be  angles  for  the  C-O-C  chain  configuration  are  in  the 
range  of  108.4°-1 16.8°.  For  the  carbon-doped  BeO  system,  the  carbon 
hybridization  bond  angles  are  very  different  from  the  standard  sp2  (120°)  and  sp3 
(109.4°)  hybridization  bond  angles,  indicating  a  higher  degree  of  bonding 
complexity.  When  one  equatorial  boron  atom  neighboring  a  C-B-C  chain  is 
replaced  by  carbon,  a  Ce-Co  bond  forms  between  the  equatorial  carbon  and  C-B-C 
chain.  The  Be-Ce-Be  bond  angle  for  this  C  is  1 10.64°,  where  as  the  Ce-Co-Be  bond 
angles  are  118.90°.  For  the  other  carbon  atom  in  the  chain  the  Be-Co-Be  bond 
angles  are  1 15. 63°-l  16.52°.  If  instead,  the  carbon  atom  replaces  a  polar  boron 
neighboring  a  C-B-C  chain  (CP-CBC),  the  Be-C-Be  bond  angles  are  in  the  range  of 
1 15.3 1°-1 16.87°  for  one  C-atom  in  the  C-B-C  chain  and  1 15. 76°-l  16.86°  for 
other  C  atom. 

To  understand  the  charge  transfer  in  the  doped  systems,  we  calculated  the  charge 
density  difference  by  subtracting  the  charge  density  of  the  isolated  C-B-C  or  C-C- 
C  chains  from  the  total  charge  density  of  the  combined  system.  The  excess  and 
depletion  charge  densities  for  this  system  indicate  that  charge  is  transferred  from 
the  B12  icosahedra  to  the  carbon  atoms  in  the  chain.  For  the  C-B-C  chain 
configuration,  the  central  boron  atom  also  transfers  charge  to  the  2  carbon  atoms 
in  the  chain.  However,  for  the  C-C-C  chain  configuration,  there  is  charge  transfer 
from  the  2  terminal  carbon  atoms  to  the  interstitial  C  atom. 

The  total  DOS  (see  Fig.  2  and  Appendix  B)  gives  further  infonnation  about  the 
nature  of  bonding  in  the  undoped  and  carbon-doped  BeO  systems.  Our  results 
indicate  that  BeO  is  a  semiconductor  with  the  calculated  band  gap  of 
approximately  2.6  eV.  The  band  gap  is  slightly  higher  than  that  of  previous 
reported  value  of  2.4  eV. 20,21  Density  functional  theory  (DFT)  is  known  to 
underestimate  the  band  gaps.  The  band  gap  reported  for  BeO  based  on  many  body 
perturbation  theory  is  3  eV22,  whereas  density  functional  tight  binding  theory  (DF- 
TBT)  calculates  a  3.45  eV  band  gap.23  For  the  pure  BeO,  the  states  near  the  higher 
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energy  part  of  the  valence  band  are  mainly  boron  2p  states,  followed  by  oxygen 
2p  and  boron  2s  states  in  the  middle  and  lower  energy  range.  The  lower  end  of  the 
conduction  band  contains  empty  boron  2p  states. 
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Fig.  2  Schematic  representation  of  total  DOS  for  pure  and  doped  B60 


Appendix  B  contains  the  total  and  atom  resolved  partial  density  of  states  for 
carbon  substituting  for  different  sites  within  the  BeO  unit  cell.  When  one  oxygen 
atom  is  replaced  by  carbon,  a  new  peak  appears  1.2  eV  above  the  Fermi  level. 
This  peak  has  a  significant  contribution  from  the  empty  2p  orbitals  of  the  carbon 
atom.  When  a  second  carbon  atom  is  doped  at  the  neighboring  oxygen,  the  DOS 
shows  significant  variations  in  the  lower  energy  part  of  the  conduction  band. 
There  are  also  some  delocalized  states  near  the  Fermi  level.  The  empty  2p  levels 
of  each  Co  atom  split  into  slightly  lower  (0.71  eV)  and  higher  energy  states  (1.56 
eV)  than  the  1.2  eV  peak  present  when  only  a  single  oxygen  is  substituted  by 
carbon.  For  the  C-B-C  chain  configuration,  no  such  peak  is  observed.  For  the  C- 
B-C  configuration,  a  new  peak  is  created  at  the  Fermi  level  with  contributions 
from  2p  orbitals  of  the  carbon  and  boron  atoms  in  the  C-B-C  chain.  For  the  C-C-C 
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chain  configuration,  the  DOS  shows  a  peak  at  1.55  eV  above  the  Fenni  level.  This 
peak  is  dominated  by  the  2p  states  of  the  interstitial  carbon  atom.  When  an 
additional  carbon  atom  is  doped  in  the  C-B-C  doped  configuration,  the  DOS  shifts 
slightly  towards  the  lower  energy  range.  This  shifting  of  DOS  towards  the  lower 
energy  is  an  indication  of  the  improved  stability  of  the  CP-CBC  and  Ce-CBC 
systems. 


3.2  Structural  and  Elastic  Properties 


The  finite  strain  method  is  the  most  common  method  to  compute  the  stiffness 
coefficients  of  materials.  In  this  approach,  the  ground  state  structure  is  strained 
according  to  symmetry  dependent  strain  patterns  with  varying  amplitudes.  The 
stress  tensor  is  detennined  after  a  reoptimization  of  the  internal  structure 
parameters.  The  elastic  stiffness  coefficients  follow  from  Hooke’s  law24,25  as  the 
proportionality  coefficients  relating  the  applied  strain  to  the  computed  stress. 


cr.  =  V6  c  £ .  . 

'  Z— </,y= 1  J  J 


(1) 


Equation  1  consists  of  a  set  of  6  linear  equations  with  27  variables:  6  components 
of  stress  and  21  elastic  constants.  The  elastic  behavior  of  a  completely 
asymmetric  material  (triclinic)  is  specified  by  21  independent  elastic  constants, 
whereas  an  isotropic  material  (cubic)  requires  only  3.  In  between  these  limits,  the 
necessary  number  is  determined  by  the  symmetry  of  the  material.  For  a  hexagonal 
crystal,  there  are  6  unique  symmetry  elements  cn,  cn,  C13,  C33,  C44,  and  C66. 
However,  since  C66=  (cn-ci2)/2,  there  are  only  5  truly  independent  symmetry 
elements.  In  Table  1,  we  compared  the  calculated  lattice  parameters  and  elastic 
constants  for  undoped  BeO  with  the  earlier  results  from  the  scientific  literature. 
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Table  1  Summary  of  the  calculated  BeO  lattice  parameters  (A)  and  single  crystal  elastic 
constants  Cij  (GPa)  with  previously  reported  values 


Property 

PAW-GGA20 

PAW-LDA26 

GGA27 

LDA28 

Our  Calculations 

lxlxl 

2x2x1 

a 

5.391 

5.328 

5.331 

5.127 

5.392 

10.792 

c 

12.312 

12.147 

12.124 

12.327 

12.312 

Cll 

603 

615 

614.99 

586.88 

585.70 

586.05 

C33 

459 

474 

474.12 

443.87 

458.99 

466.06 

Cl2 

109 

140 

122.22 

133.19 

128.18 

121.66 

Cl3 

50 

61 

46.51 

47.41 

53.30 

45.29 

C44 

179 

180 

192.38 

176.87 

178.38 

C14 

22.75 

22.4 

22.25 

The  bulk  modulus  B  and  shear  modulus  G  can  be  measured  to  detennine  elastic 
properties  of  materials.  On  the  basis  of  the  scheme  of  Voigt29  (assuming  uniform 
strain)  and  Reuss,30  (assuming  unifonn  stress)  and  with  the  approximation  of 
Hill’s31  arithmetic  average,  we  calculated  the  bulk  modulus  B  and  shear  modulus 
G  for  B60.  Tables  2,  3,  and  4  compare  the  calculated  lattice  parameters,  elastic 
constants,  and  elastic  moduli  for  the  pure  and  doped  BeO  systems.  Furthermore, 
the  Young  modulus  E  and  Poisson’s  ratio  a  have  been  calculated  from  the 
following  equations  for  the  isotropic  material: 

E_  9 BG  _  (3B-2G) 

~  (3B  +  G)’a  ~  2(3B  +  G)' 


Table  2  Calculated  lattice  parameters  (A),  cell  volume  (A3)  and  single  crystal  elastic 
constants  Cij  (GPa)  for  pure  and  carbon-doped  BsO 


Property 

b6o 

C(®,BP 

C(a)0 

2C(a)0 

C-B-C 

c-c-c 

C-B-C  +  CP 

a 

10.792 

10.791 

10.806 

10.819 

10.833 

10.832 

10.825 

c 

12.312 

12.307 

12.307 

12.307 

10.833 

12.302 

12.297 

V 

1,241.94 

1,241.29 

1,244.56 

1,247.65 

1,250.82 

1,250.04 

1,248.28 

Cll 

586.05 

578.21 

582.24 

583.07 

581.6 

579.1 

585.1 

C33 

466.06 

459.2 

460.89 

466.73 

472.3 

469.6 

472.1 

C12 

121.66 

121.68 

117.52 

119.38 

120.1 

120.8 

122.1 

Cl3 

45.29 

44.54 

42.04 

48.94 

51.6 

49.4 

51.4 

C44 

178.38 

177.92 

176.43 

172.56 

171.6 

173.5 

176.0 

C14 

22.25 

20.25 

19.87 

19.35 

15.9 

21.0 

21.6 
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Table  3  Calculated  lattice  parameters  and  elastic  constants  Cij  (GPa)  for  BeO  with  2  C 
atoms  substituting  for  neighboring  O,  separated  by  O  interstitial  (Co-Oi-Co  chain) 


a 

b 

c 

Cn 

C12 

C13 

C14 

C15 

10.810 

10.834 

12.316 

540.88 

79.98 

8.32 

3.94 

5.78 

Cl6 

C22 

C23 

C24 

C25 

C26 

C33 

C34 

4.27 

530.25 

5.14 

4.84 

24.68 

5.58 

410.57 

5.7 

C35 

C36 

C44 

C45 

C46 

Css 

C56 

C66 

7.30 

5.78 

210.57 

0.89 

44.93 

180.38 

0.96 

126.38 

Table  4  Comparison  of  the  calculated  elastic  moduli  (V-Voight,  R-Reuss,  H-Hill)  in  GPa, 

Poison’s  ratio,  and  polycrystalline  elastic  anisotropy  Voight  Reuss  (AYR) 

Property 

BsO 

C(®BP 

c®o 

2C(«),0 

Co-Oi-Co 

C-B-C 

c-c-c 

C-B-C  +  CP 

By 

229.19 

226.35 

225.40 

229.71 

185.39 

231.38 

229.65 

232.45 

Gv 

212.84 

210.48 

211.96 

209.77 

196.01 

208.93 

209.13 

211.2 

br 

223.31 

220.41 

219.47 

224.33 

178.12 

226.5 

224.6 

227.27 

Gr 

205.86 

204.2 

205.11 

202.64 

177.94 

202.35 

202.12 

204.28 

Bh 

226.25 

223.38 

222.43 

227.02 

181.76 

228.94 

227.12 

229.86 

Gh 

209.35 

207.34 

208.54 

206.21 

186.98 

205.64 

205.62 

207.74 

Eh 

480.01 

475.05 

476.65 

474.85 

417.7 

474.77 

473.87 

478.94 

aH 

0.15 

0.15 

0.14 

0.15 

0.12 

0.15 

0.15 

0.15 

AYR  (%) 

1.668 

1.513 

1.642 

1.727 

4.831 

1.601 

1.704 

1.666 

To  the  best  of  our  knowledge,  there  are  few  experimental  studies  of  the  elastic 
constants  for  BeO.  Petrak  et  al.32  used  the  resonant-sphere  technique  to  obtain  the 
aggregate  elastic  moduli  for  hot-pressed  BeO  specimens  under  ambient 
conditions.  He  reported  the  bulk  modulus  B,  shear  modulus  G,  and  Young’s 
modulus  E  as  230,  206,  and  472  GPa,  respectively,  whereas  Andrievskiy33 
reported  228,  204,  and  471.4  GPa,  respectively.  Earlier  ab  initio  studies  report 
that  the  bulk  moduli  (B)  for  BeO  vary  from  222  to  235  GPa.34-37  Our  calculated 
value  of  Bh  =  226.25  GPa  agrees  with  the  previous  studies.  Our  calculated  bulk 
modulus,  shear  modulus,  and  Young’s  modulus  are  also  in  agreement  with  the 
experimental  findings.33,34 

Carbon  doping  significantly  changes  both  the  elastic  constants  and  elastic  moduli. 
The  doping  of  one  carbon  atom  at  a  boron  or  oxygen  site  decreases  the  elastic 
moduli  between  0.4%  and  1.7%.  The  doping  of  2  carbon  atoms  at  2  oxygen  sites 
increases  the  bulk  modulus  by  0.3%,  yet  decreases  the  shear  and  Young’s 
modulus  by  0.5%.  When  there  is  an  additional  interstitial  oxygen  between  the  2 
Co,  the  elastic  moduli  are  significantly  decreased.  For  this  case,  the  bulk  modulus 
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decreases  by  20%,  whereas  the  shear  and  Young’s  moduli  decrease  by  10%— 13%. 
It  is  interesting  to  note  the  change  in  the  elastic  moduli  for  C-B-C  and  C-C-C 
local  chain  configurations  within  B60.  When  either  a  boron  or  carbon  atom  is 
doped  as  the  interstitial  between  2  Co,  the  bulk  modulus  increases,  whereas  the 
shear  modulus  and  Young’s  moduli  decrease.  The  C-B-C  chain  configuration  has 
higher  B,  G,  and  Y  moduli  than  the  C-C-C  chain  configuration.  If  an  additional 
carbon  is  doped  in  either  a  polar  or  equatorial  sites  of  a  icosahedra  neighboring  a 
C-B-C  chain,  further  improvements  in  the  3  moduli  can  be  seen.  However,  the  G 
and  Y  values  are  lower  than  that  of  the  undoped  BeO. 

4.  Conclusions 


In  summary,  we  calculated  the  structural  and  elastic  properties  of  pure  and 
carbon-doped  Bf,0  under  ambient  conditions  within  the  plane-wave 
pseudopotential  DFT.  We  demonstrate  how  carbon  doping  can  be  used  to  modify 
the  electronic  structure  and  elastic  properties  of  the  hexagonal  B60.  It  is 
energetically  preferable  for  carbon  to  replace  oxygen  rather  than  boron.  When 
carbon  replaces  oxygen,  the  lattice  parameters  and  cell  volume  expand  to  relieve 
the  residual  stress.  The  energetic  cost  of  carbon  doping  increases  when  2  carbon 
atoms  replace  2  oxygen  atoms.  Adding  an  interstitial  oxygen  along  a  Co-Co  chain 
significantly  changes  the  B60  unit  cell  and  decreases  the  elastic  moduli  leading  to 
a  reduction  in  material  hardness.  In  contrast,  interstitial  boron  or  carbon  doping 
leading  to  local  C-B-C  or  C-C-C  chain  configurations  improves  the  bulk  modulus 
at  the  expense  of  the  shear  and  Young’s  modulus.  The  C-B-C  chain  configuration 
results  in  higher  B,  G,  and  Y  moduli  than  the  C-C-C  chain  configuration.  The  3 
moduli  further  improve  as  additional  carbon  atoms  are  substituted  for  boron  in 
either  the  polar  or  equatorial  site  of  a  neighboring  B12  icosahedra.  These 
calculations  suggest  that  doping  can  either  harden  or  soften  BeO.  The  hardening 
or  softening  depends  on  the  location  and  the  type  of  a  dopant.  When  the  doping 
creates  a  B/jC-like  local  structure  within  B60,  there  is  an  increase  in  the  bulk 
modulus  and  decrease  in  shear  and  Young’s  moduli.  This  study  demonstrates  that 
the  B/tC-like  local  doping  could  be  used  to  improve  the  mechanical  stability  and 
elastic  properties  of  BeO. 
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Appendix  A.  Electron  Localization  Function  for  Pristine  and 

Carbon-Doped  BeO 
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Fig.  A-l  Two-dimensional  (2-D)  plots  for  total  charge  density  (a-d)  and  electron  localization 
function  (ELF)  (e-h)  for  the  B60.  The  plane  normal  to  |001]  direction:  a-c  and  e-g;  the 
plane  normal  to  |110]  direction:  d  and  h. 
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(a)  1C@0 


(b)  2C@0 


(d)  C-O-C 


(e)  C-C-B 


(f)  C-B-C 


Fig.  A-2  2-D  plots  for  the  total  charge  density  and  ELF  for  C-doped  BeO  case:  a)  the  plane 
normal  to  |001]  direction;  b)  and  f)  the  plane  normal  to  1110]  direction. 
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Charge  Density 


ELF 


(c)  C-B-C  +  CE 


(d)  C-B-C  +  CE 

Fig.  A-3  2-D  plots  for  the  total  charge  density  and  ELF  for  C-  doped  BsO  case,  a)  represents 
the  plane  containing  doped  C-B-C  atoms  and  Cp  atom;  b)  represents  the  plane  normal  to 
[110]  direction;  3c)  represents  the  plane  containing  doped  C-B-C  atoms  and  Ce  atom;  3d 
represents  the  plane  normal  to  [110]  direction. 
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Appendix  B.  Atoms  Resolved  Density  of  States  for  Pristine  and 

Carbon-Doped  BeO 
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Fig.  B-l  Total  density  of  states  (TDOS)  and  atom  resolved  (O,  Be,  and  BP)  partial  density  of 
states  (PDOS)  for  B60 
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Fig.  B-2  TDOS  and  atom  resolved  (Co,  nearest  Be)  PDOS  for  BsO  with  one  C-atom 
substituted  for  O 
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Fig.  B-3  TDOS  and  atom  resolved  (Co,  nearest  Bc)  PDOS  for  660  with  2  C  substituted  for 
neighboring  O.  Both  Co  have  similar  PDOS. 
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Fig.  B-4  TDOS  and  atom  resolved  (Co,  Oi,  nearest  Bc)  PDOS  for  C-O-C  chain  substitution 
within  BfiO.  Both  Co  have  similar  PDOS. 
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Fig.  B-5  TDOS  and  atom  resolved  (Co,  Ci,  nearest  Be)  PDOS  for  C-C-C  chain  substitution 
within  B6Q.  Both  Co  have  similar  PDOS. 
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Fig.  B-6  TDOS  and  atom  resolved  (Co,  Bi,  nearest  Be)  PDOS  for  C-B-C  chain  substitution 
within  B6Q.  Both  Co  have  similar  PDOS. 
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Fig.  B-7  TDOS  and  atom  resolved  (Co,  Bo,  Ci)  PDOS  for  C-C-B  chain  substitution  within 

b6o 
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Fig.  B-8  TDOS  and  atom  resolved  (Co,  Bi,  Ce)  PDOS  for  Ce  +  C-B-C  chain  substitution 
within  BfiO.  Cl  is  bonded  to  Ce,  C2  is  bonded  to  Be. 
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Fig.  B-9  TDOS  and  atom  resolved  (Co,  Bi,  Cp)  PDOS  for  Cp,  +  C-B-C  chain  substitution 
within  B6Q.  Cl  is  separated  from  CP  by  lBe,  C2  is  bonded  to  Be. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


AVR 

anisotropy  Voight  Reuss 

DOS 

density  of  states 

DF-TBT 

density  functional  tight  binding  theory 

DFT 

Density  functional  theory 

ELF 

electron  localization  function 

GGA 

generalized  gradient  approximation 

TDOS 

total  density  of  states 

PAW 

projector  augmented  wave 

PDOS 

partial  density  of  states 

VASP 

Vienna  Ab  initio  Simulation  Package 

29 


1  DEFENSE  TECHNICAL 
(PDF)  INFORMATION  CTR 

DTIC  OCA 

2  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIO  LL 

IMAL  HRA  MAIL  &  RECORDS 
MGMT 

1  GOVT  PRINTG  OFC 
(PDF)  AMALHOTRA 

1  DIRUSARL 
(PDF)  RDRL  WMM  E 
J  S  DUNN 


30 


